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1. Introduction 

The starting formula in calculations of the Casimir energy Eq = (1/2) ^u;„ requires the 
knowledge of the spectrum of the quantum field model under consideration [T] . For this 
aim the boundary conditions at the spatial infinity should be imposed on the field in the 
case of an unbounded configuration space. Making here use of the radiation condition 
leads, as a rule, to a complex discrete energy spectrum (open systems with energy losses 
due to the outgoing waves). In this situation a more admissible is the condition utilized 
in formulation of the scattering problem, i.e., at the spatial infinity the filed should be 
a sum of incoming and outgoing waves. The spectrum of the same Hamiltonian with 
scattering conditions is nonnegative but continuous, and the respective natural modes 
are not squared integrable. In order to define a correct integration over such a spectrum 
one has to introduce the spectral density, which is expressed in terms of the Jost function 
of the relevant scattering problem. The frequency equation, which is derived by making 
use of the radiation condition, is related, in a direct way, with the Jost function. As 
a result the final formulas for the spectral sums and spectral integrals turn out to be 
identical in the end. Thus it is proposed a rigorous justification of the results of Casimir 
calculations obtained by employment of the radiation condition and frequency equations 
with discrete set of complex roots. 

2. Complex frequencies and quasi-normal modes in unbounded oscillating 
regions 

Here we show in the general case in what way complex frequencies and quasi-normal 
modes appear when considering harmonic oscillations in unbounded regions. Let a 
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closed smooth surface S divides the d-dimensional Euchdean space M"^ into a compact 
internal region Di^ and noncompact external region D^^. We consider here a simple 
scalar wave equation 
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A-^^J«(t,x) = 0, (1) 

where c is the velocity of oscillation propagation and A is the Laplace operator in M°'. 
For harmonic oscillations 

M(t,x) = e-'"*M(x) (2) 

the wave equation ((T)) is reduced to the Helmholtz equation 

(A + m(x) = 0, A; = ^ . (3) 

The oscillations in the internal region Din are described by an infinite countable set 
of normal modes 

u„(t,x) =e-^""V(x), n=l,2,.... (4) 

The spatial form of the normal modes (the functions m„(x)) is determined by the 
boundary conditions which are imposed upon the function m(x) on the internal side 
of the surface 5*. These conditions should fit the physical content of the problem under 
study. The set of normal modes is a complete one. Hence any solution of © obeying 
relevant boundary conditions can be expanded in terms of the normal modes Mn(x). 

When considering the oscillations in the external domain Dcx one imposes, in 
addition to the conditions on the compact surface S, a special requirement concerning 
the behavior of the function m(x) at large r = |x|. Usually for this purpose the radiation 
boundary conditions, proposed by Sommerfeld [21 E] , are used 

lim r^~u{r) = const , lim ( — iku \ = . (5) 



For real values of the wave vector k (for real frequencies u) the solution to 0, which 
obeys the radiation conditions © and reasonable boundary condition on a compact 
surface S, identically vanishes. 

If we remove the requirement of reality of the wave vector k, then the statement, 
formulated above, does not hold. Namely, the wave equation and the Helmholtz 
equation Q have nonzero solutions with complex frequencies, these solutions obeying 
the radiation boundary conditions ^ and a common condition on a closed compact 
surface S (for instance, Dirichlet or Neumann conditions). 

As a very simple and clear example we consider the oscillations of electromagnetic 
field outside a perfectly conducting sphere of radius a. In this case the electric and 
magnetic fields are expressed in terms of two scalar functions fki^{r) and f™{r) (Debye 
potentials) which are the radial parts of the solutions to the scalar wave equation (^. 
Outside the perfectly conducting sphere placed in vacuum the solution to the Helmholtz 
equation ^ obeying the radiation conditions (0) has the form {d = 3) 

/fcKO = C/i;'^(^r), r>a, (6) 
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where hf^^{z) is the spherical Hankel function of the first kind At the surface of 
perfectly conducting sphere the tangential component of the electric field should vanish. 
This leads to the following frequency equation for TE-modes 

/^«(^^a)=0, />1 (7) 

and for TM-modes 

±(fti"(^r))=0. r = a. />1. (8) 

The spherical Hankel function hf^^{z) is e^^ multiplied by the polynomial in 1/z of a 
finite order. Hence frequency equations (|7j) and (jH)) have a finite number of roots, which 
are in the general case complex numbers. For / = 1 (the lowest oscillations) equations 
dlj) and (jHI) assume the form {z = au/c) 

d 



h\'>{z) =--e'Ml + -)=0 (TE modes), (9) 



— (^z hf~^ {z)^ = - ^ e^' {z^ + iz - 1) = (TM modes) . (10) 
Thus the lowest eigenfrequencies are 

- = - - (TE modes), (11) 
c a 

-= - — (i ± V3) (TM modes) . (12) 
c 2a 

The complex eigenfrequencies lead to a specific time and spatial dependence of the 
respective natural modes and ultimately of the electromagnetic fields. So, with allowance 
of (fTTjl . we obtain 

e-'"VfcT(r) = -iC-e('-^*)/'^fl--V r>a. (13) 

Thus, the eigenfunctions are exponentially going down in time and exponentially 
going up when r increases. Such time and spatial behaviour is typical for eigenfunctions 
describing oscillations in external unbounded regions, the physical content and details 
of oscillation process being irrelevant. The eigenfunctions corresponding to complex 
eigenvalues are called quasi-normal modes keeping in mind their unusual properties 0. 
The physical origin of such features is obvious, in fact we are dealing here with open 
systems in which the energy can be radiated to infinity. Therefore in open systems field 
cannot be in stationary state. 

Quasi-normal modes do not obey the standard completeness condition and the 
notion of norm cannot be defined for them Therefore these eigenfunctions cannot 
be used for expansion of the classical field with the aim to quantize it and to introduce 
the relevant Fock operators. Thus the radiation conditions are not appropriate for the 
Casimir calculations. 

The same situation, with regard to quasi-normal modes, takes place when we 
consider the oscillations of compound media. In this case in both the regions D^^ and 
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Dcx the wave equations are defined 
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with the matching conditions at the interface S, for example, of the following kind 



where nin(x) and nex(x) are normals to the surface S at the point x for the regions 
An and Ax 5 respectively. The parameters Cin, Cgx, Kn, and Aex specify the material 
characteristics of the media. At the spatial infinity the solution Mex(t,x) should satisfy 
the radiation conditions ©. For real k we again have only zero solution in this problem, 
both functions Min(t, x) and Mcx(i,x) vanishing. However, the wave equations (fT^ and 
(|15p have nonzero solutions with complex frequencies, i.e. quasi-normal modes, which 
satisfy the matching conditions (fTBj) and (fTTj) at the interface S and radiation conditions 
(0) at spatial infinity. It is important, that the frequencies of oscillations in internal 
(Din) and external (-Dex) regions are the same. A typical example here is the complex 
eigenfrequencies of a dielectric ball. This problem has been investigated by Debye in his 
PhD thesis concerned with the light pressure on a material ball jH|. 

The physical content of the radiation conditions is very clear. They select only the 
oscillations with real frequencies caused by external sources, which are situated in a 
compact spatial area. From the mathematical point of view, these conditions ensure the 
uniqueness of the solution of the nonhomogeneous boundary problems in the external 
region D^x or in the whole space (Din + -Dcx) in the case of compound media. 

When formulating the radiation conditions in the text books, only the real wave 
vector k is considered. The possibility of existence of quasi-normal modes with complex 
frequencies satisfying the radiation conditions at spatial infinity with complex wave 
vector k is not mentioned usually. I am aware only of one literature source where 
the eigenf unctions with complex frequencies are noted in this context. It is the article 
written by Sommerfeld in the book [Hj. 

The standard physical method to escape appearance of complex frequencies and 
quasi-normal modes is the following: the system under study is placed in a sphere of 
a large radius R and the field is subjected to an additional boundary condition on 
this sphere. The initial differential operator has now a real discrete spectrum, and the 
respective eigenfunctions are normalized in L2. Upon conducting the calculations the 
limit R 00 should be taken. 

However, in order to justify this approach one has to prove that the final result does 
not depend on the explicit form of auxiliary boundary conditions imposed on the large 
sphere r = R. Obviously, it is nontrivial task. Furthermore, in this approach always 
the feeling is rest, that something is lost because the region r > R is actually ignored. 




(16) 



(17) 
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Fortuitous, there is a rigorous way to consider the open systems by making use of the 
formahsm of the scattering problem without introducing a large auxiliary sphere. 

3. Description of open systems by scattering method 

In the scattering problem the same differential operator (for example, the Laplace 
operator Q) is considered, but instead of the radiation conditions (0) it is required 
that the eigenf unctions (f)i{k,r) at spatial infinity r ^ oo are reduced to the sum of 
incoming and outgoing waves 

Mk,r)^-^[Fii-k)e-"''- -Fi{k)e"''-], r oo . (18) 

Here Fi{k) is the Jost function, and spherical symmetry in the problem under 
consideration is expected. The operators of the type 

-A + V{r) (19) 

have real continuous spectrum 

0<P <oo, (20) 

when the potential V{r) obeys known conditions at r = and at infinity r —>■ oo. 

Due to the continuity of the spectrum the corresponding eigenfunctions are not 
normalized in L2. Making here use of the normalization on the Dirack ^-function proves 
to be useful in many physical applications, however it does not allow one to define, in 
a correct way, the integration over a continuous spectrum. Indeed, let ipy.i'x) be the 
eigenfunctions corresponding to the spectrum (pUj) normalized on the Dirack (5-function 

j V'k(x)V^k'(x) d^x = 5(3) (k - k') , (21) 

for example, those can be plane waves 

*W=(^»*" (22) 
Due to its meaning, the spectral density p{k) should be defined in the following way 

p{k) = 1 4^1{^)M^) d^x = 5(3) (0) = (23) 

where V is the volume of the space and obviously ^ 00. It is worth noting 
that the spectral density (j2SI), considered just on the formal footing, does not depend 
on the details of the dynamics, i.e. on the potential V{r) and corresponding boundary 
conditions, therefore it can not be used in physical calculations. 

In the rigorous mathematical theory of scattering processes this problem has been 
solved. In the most simple way the spectral density for scattering states can be derived 
as a consequence of a spectral property j7]. Let r{z) = {h — z)~^ and tq^z) = {ho — z)~^ 
be the resolvents of a complete (h) and free [ho) Hamiltonians. When the potential 
meets the relevant conditions, then the relation, known as the spectral property, holds 

2 Im Tr [r{E + iO) - ro{E + iO)] = tr [q{E)] , (24) 



Boundary conditions at spatial infinity 



6 



where 



q{E) = -iS\E)^S{E) (25) 
is the two-body time delay operator and S is the scattering matrix 

5(A:) = e-^« = ^. (26) 

Here 6{k) is the phase shift and F{k) is the Jost function. 

By making use of the spectral representation for the resolvents and definition (j^^ 
one obtains immediately 

Ap{k)^p{k)-po{k) = ^^^\nS{k) 

It is worth noting that in this formula the contribution of a free unbounded space is 
already subtracted. The finite expression for Ap{k) is a direct consequence of considering 
complex values of the wave vector k and respective energy k'^. 

4. Correct transition to imaginary frequencies in the case of unbounded 
domains 

In the Casimir calculations the spectrum of quantum field systems is not known 
explicitly. As a rule, it is given by the roots (un) of a frequency equation 

fico) = (28) 

with the known function /. In order to sum these roots the counter integration can be 
used 

where the counter C encloses, counterclockwise, all the roots of 

For a compact configuration space without dissipation the roots of ()28p are real and 
positive. The counter C can be deformed continuously into the counter Cr with R oo 
(see figure la). Usually f{x) = f{Vx^), therefore f{x) = f{—x). The integral along the 
semicircle of radius R vanishes when i? ^ oo or it is canceled by the subtraction of the 
free space contribution. As a result we obtain 

Let us derive the analog of (jHUjl for the open systems. According to the inferences 
of the previous section we have to consider the following integral along the continuous 
spectrum 
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where F{k) is the Jost function in the pertinent scattering problem. This function is 
analytic in the lower half plane of the complex variable k. We assume that there are 
no bound states in the system in question. Under this assumption the zeros of the Jost 
function F{k) can lay only in the upper half plane IHUHl- In view of this we infer that 
the counter integrals (see figure lb) 

vanish because inside the counter (C/j) there are no zeros and singularities of the 
function F{—k) {F{k)). Now we can transform in ()31|) the integration along the real 
positive axes k to the integration along the imaginary frequencies k = iy, y > with 
the result 

- kAp{k)dk = — I lnF(-iy)dy. (33) 
2 Jo TT Jo 2 dy 




Figure 1. The contours on the complex k plane which are used when going on to the 
integration over the imaginary frequencies in the case of bounded (a) and unbounded 
(b) configuration spaces. 



Now we take into account the following. The function f{uj) in frequency equation 
(I28|l ■ which is derived by making use of the radiation conditions at spatial infinity, is 
actually the Jost function up to irrelevant multiplier. Indeed, the solution (fTHj) will 
contain only outgoing waves when 

F{-k)=0. (34) 

For simplicity we have dropped here the index /. Thus we have 

f{uj) = F{-k), e = u?lc^. (35) 

Now we can infer that the right-hand sides of ()3()|1 and coincide. Hence, the formula 
(IHn|) with the function / defining the frequency equation is applicable also to unbounded 
configuration spaces. 
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All this removes the objections jTUl E] against using for unbounded regions 
and explains why it gives correct results there pp. 

5. Conclusion 

Importance of proper choosing conditions, imposed on the field at spatial infinity when 
determining the physical spectrum in Casimir calculations, is illustrated by a circular 
infinite dielectric cylinder. If we confine our consideration to the oscillating asymptotics 
for r ^ oo the spectrum of Maxwell equations is continuous. Only this situation has 
been considered in all pertinent papers (see, for example, ^3). However from the 
physical point of view are acceptable also solutions exponentially going down when 
r — > oo. It is these configurations of electromagnetic field (surface waves with discrete 
spectrum) that are used in dielectric wavegaudes [iH] • Thus the question arises whether 
this branch of spectrum has been taken into account in relevant calculations. Casimir 
studies, which allow for the material properties of the boundaries, show importance 
of the surface waves, or more precisely, of surface plasmon contributions ^1] at short 
distances (regime without retardation). Obviously, such bound states should be added 
to the spectral density 
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